Charge renormalization for effective interactions of colloids at water interfaces 
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We analyze theoretically the electrostatic interaction of surface-charged colloids at water interfaces 
with special attention to the experimentally relevant case of large charge densities on the colloid- 
water interface. Whereas linear theory predicts an effective dipole potential the strength of which is 
proportional to the square of the product of charge density and screening length, nonlinear charge 
renormalization effects change this dependence to a weakly logarithmic one. These results appear 
to be particularly relevant for structure formation at air-water interfaces with arbitrarily shaped 
colloids. 
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The effective interactions of colloids trapped at fluid 
interfaces reveal qualitatively new features when com- 
pared to the ones in colloidal bulk solutions. First, there 
is the possibility of long-ranged capillary attractions me- 
diated by deformations of the interface [l|, 0|. Second, 
many colloids carry a significant amount of charge (e.g. 
charge-stabilized polymeric colloids, mineralic disks, pro- 
teins) and the exponentially screened electrostatic inter- 
actions in ionic bulk solvents become longer-ranged at in- 
terfaces between water and a nonpolar medium (typically 
air or oil) . At such interfaces the colloids exhibit effective 
dipole-like repulsions which lead to the stabilization of 
two-dimensional crystals even at low surface coverages 
These effective dipoles originate from colloidal sur- 
face charges on the water side and a cloud of screening 
ions in the water phase which is asymmetric with respect 
to the interface plane. Within a simple model (the only 
analytically tractable one) the colloids are approximated 
as equal point charges q located in the interface plane 
and the water phase is treated as a linearly screening 
medium. To leading order the interaction between two 
chafes q in the interface plane at separation d is given 

U{d)^q'-^^^. (1) 

Here, ei and e2 are the permittivities of the nonpolar 
medium and water, respectively, and eq is the dielectric 
constant of vacuum. According to this linear model the 
repulsion depends quadratically on the Debye screening 
length = (e2eo/(2/3coe^))^/'^ where cq is the concen- 
tration of monovalent ions in bulk water, e is the elemen- 
tary charge, and /3^^ = fc^T. On this basis one would 
expect the repulsion U oc c^^ to become significantly 
weaker upon adding electrolytes. Various studies of col- 
loidal aggregation at interfaces have used the predictions 
of the linear model for quantitative analysis of experi- 



mental results (see, e.g., Refs. @, jij). In Ref. [1| g in 
Eq. Il]) was replaced by q^s oc to account for geo- 
metric effects of a charged spherical colloid, which leads 
to J7 oc c^^. 

The high colloidal surface charge densities CTc on the 
water side of experimentally used colloids (easily up to 
0.5 e/nm^) invalidate the naive use of the linearized 
Debye-Hiickel (DH) model with bare charges. Strong 
charge renormalization will occur due to the nonlinear 
contributions of the governing Poisson-Boltzmann equa- 
tion (PB) in the water phase. The renormalization pro- 
cedure (based on the separation of length scales) consists 
of the identification of the appropriate corresponding lin- 
ear solution of the PB problem at distances > from 
the charges. There the electrostatic potential $ is small 
and linear DH electrostatics holds: ~ For a 

uniformly charged wall or sphere, this solution has the 
same functional form as if the entire problem is solved 
within the linear theory and the nonlinear effects alter 
only the prefactor. This prefactor leads to a renormal- 
ized, effective charge Q. For non-spherical charged bod- 
ies the map between the DH solution and the PB solution 
in the linear region requires a selection of the appropriate 
boundary conditions at the charged object such that the 
DH and the PB solution match at the far field Q . In the 
limit (Tc ^ oo of the surface charge density the renor- 
malized DH potential at the colloid surface levels off at a 
constant regardless of the geometry of the charged body. 

The renormalization of charges at an interface is ex- 
pected to differ from that in the bulk due to the proxim- 
ity of a nonpolar phase which induces an algebraic decay 
of the electrostatic field near the interface; to a large ex- 
tent its strength is determined by the potential within the 
screening length. In order to study the effect of an inter- 
face on the renormalization we have chosen the experi- 
mentally relevant system of a charge-stabilized colloidal 
sphere trapped at an interface with water. The renor- 
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<P(z-~-oo) ~ exp(KZ) / z 



FIG. 1: Side view of a single colloid (homogeneously charged 
on the water side) trapped at the interface. Most of the 
counterions are confined in a layer close to the colloid sur- 
face with a width of the order of the Gouy-Chapman length 
~ 2e2eo/(/3e(Jc). In many colloidal experiments, Zg(~ 1 
nm) < K'^i^ 1...300 nm) < 1 pm). 



malized dipole field can be described in terms of a single 
renormalized parameter given by the effective charge goff- 
We find that the ratio qcs/q factorizes into a geometric 
part (describable by a linear theory) which takes into ac- 
count the geometry of a charged object, i.e., the charge 
distribution at the colloid-water interface, and a nonlin- 
ear part which is described by the analytically solvable 
case (within PB theory) of a charged wall, thus being in- 
dependent of the colloid shape and the contact angle. In 
contrast to the case of colloids in the bulk we find that q^s 
does not level off for highly charged particles. Also the 
functional dependence of qcs on k differs from the bulk 
case; nonetheless geff remains an increasing function of 
K As a consequence, the effective repulsion given by 
Eq. ([1]) becomes only weakly dependent on the screening 
length. 

The model. For a single spherical colloid of ra- 
dius R trapped at an interface as indicated in Fig. [T] 
we have solved the electrostatic problem given by the 
Poisson-Boltzmann equation in the water phase, V*^0 = 
K*^ sinh[(/)], and the Laplace equation in the oil phase, 
V*^(/> = 0. Here, <p = e^$, V* = i?V and k* = kR are 
the dimensionless electrostatic potential, gradient oper- 
ator and screening length, respectively. At the water- 
oil and the colloid-oil interface the tangential electric 
field and the normal electric displacement are continuous, 
while at the colloid-water interface the normal electric 
displacement has a jump CTc. The differential equations 
with the appropriate boundary conditions are solved us- 
ing the finite element method package FEMLAB ^ld\ . In 
order to determine the potential at large distances from 
the particle we have chosen the computational space to 
be 8000 R so that the boundary conditions enclosing the 
box do not infiuence the data of interest. The nomi- 
nal charge on the colloid is g = Cc 2ttR^{1 + cos 9). 
We have determined the effective charge through equat- 
ing the asymptotics of the potential in the water-oil in- 
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FIG. 2: The renormalization function in the linear regime. 



terfacial plane to the asymptotics of the potential for 
the point charge in Debye-Hiickel approximation: q^ff = 
linip^oo(27reoe2/ei)(p^/ e/3K~^) 0(p, z — 0). Finite box 
size effects become visible at a distance p — 500 R from 
the colloid; thus all our data are taken within this range. 
The electrostatic interaction between two colloids at sep- 
aration d is indeed given by Eq. ([1]) to leading order in d, 
with q replaced by goff ■ This can be shown by a direct cal- 
culation of the force via a pressure tensor integration over 
the midplane (symmetry plane) between the two colloids. 
For fixed permitivities, the ratio qes/q = 9{k* ,<7*;0) de- 
fines a renormalization function which depends on k*, 
the dimensionless charge density a* — (e/3i?/(eoe2)) 
and 9. 

The linear Debye-Hiickel regime. The linear regime 
holds if (/> <C 1 everywhere and corresponds to a*K*~^ <^ 
1. (The retrieval of the linear regime in this limit can be 
confirmed from the exact solution for the charged wall 
model.) In this regime, the renormalization function is 
independent of cr*: g ^ 5iin('«*, 0)- The variation of gun 
with K* and 9 is moderate and thus the renormalization 
function is of the order 1 (see Fig. [5]). The variation 
of 5iin resembles a weak effective power-law for a lim- 
ited range of k* but it is clearly inconsistent with the 
proposal in Ref. [ill that it should vary oc k*~^ in the 
range 1 < k*^^ < oo. The weak dependence of gun on 
K* reflects the fact that the electrostatic field originating 
from the surface charges "escapes" to the insulator phase 
both through the colloid and, to some extent, through the 
electrolyte. At large k* the electrolyte "escape" route is 
blocked due to the thick counterion cloud surrounding the 
charged colloid and so the dependence of gun on k* disap- 
pears. The inadequate assumption of Ref. ^6*1 is that the 
electrolyte "escape" route is the only one except for the 
field originating from the charges near the three-phase 
contact line. 

The nonlinear regime. As inferred from the linear 
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FIG. 3: The charge renormalization function in the nonlinear 
regime. For a colloid of radius i? = 1 ptm, the two dimension- 
less charge densities cr* = 500 and 8000 correspond to charge 
densities of 0.9 and 15 /xC/cm^ which approximately bracket 
the charge densities occurring on polymeric colloids. 

regime the geometric contributions to the effective charge 
do not have a strong influence on g (giin is of the or- 
der of 1 for various contact angles as shown in Fig. [2]). 
This encourages us to deduce some general properties of 
g without solving the full problem explicitly. In typi- 
cal colloidal experiments 5] the radius of the colloid is 
of the order of 1 /xm and thus is much larger than the 
screening length for electrolyte concentrations cq > 10~^ 
M [kT^ < 0.1 fim). Therefore close to the colloid surface 
at the water side the electrostatic problem is similar to 
that for a charged wall in electrolyte. 

Since for a charged wall the potential outside the 
screening length levels off at large CTc 0] and the strength 
of the potential in the linear regime is g*k*~^, this im- 
plies that for fixed and large cr*, limo-*^ooff 
(in order to satisfy cff^*^^ = const), and that for 
fixed, large cr*, g must increase with k*, i.e., g must in- 
crease upon adding electrolyte. Extrapolating these re- 
sults for the charged wall to the present situation we find 
g « 4/(ct*«;*~^), i.e., geff is proportional to the screening 
length. Thus the interaction potential between two col- 
loids (Eq. ([T]) with q (/off) is independent of the screen- 
ing length and thus of the electrolyte concentration, at 
least within this crude "wall approximation" . 

Our numerical results show, however, that, unlike in 
the bulk case, the effective charge does not level off but 
increases slowly: q^s oc In cr*. This can be understood 
in terms of a second, somewhat more refined "wall ap- 
proximation". At the interface, the asymptotic behav- 
ior of the potential is determined by the electric field 
which "escapes" to the oil phase. The escaping field 
strength is proportional to the potential right at the 
colloid surface on the water side because the escaping 



field lines originate there. Thereby we can approximate 
the charge renormalization function from the contact po- 
tential at the wall gwaii = <^ces/'^*- The relation be- 
tween the surface charge and the potential at contact 
'/'wall for a charged wall is a* = 2k* sinh[0^^ji/2] 0]; 
in the linear limit (i.e., small a*/K*) this reduces to 
'?^waii ~ '^c'^*~^ and in the highly nonlinear limit (i.e., 
cr*/K* large) (|)'^^^^ = 2 ln(cr*K*~^). The effective surface 
charge is obtained by equating the two limiting cases 
leading to gwaii = 2(ln(cr*K*~^))/(CT*K*~"^). However, 
the full renormalization function g contains in addition 
the geometric contributions unaccounted for by the wall 
approximation. We augment the nonlinear "wall" part 
by the linear "geometry" part, which we have shown in 
Fig- El 5 ~ 5wan(cr*,K*)giin(K*,6')- In the strongly non- 
linear regime this ansatz describes our full numerical data 
for g rather well (see Fig. [3]). The wall model approxima- 
tion of the renormalization function can be corroborated 
in an alternative, more involved determination of cr* by 
calculating the effective dipole generated by the surface 
charges and the counterion cloud. The latter approach 
gives rise to corrections 0{k*~^) which explain the be- 
havior of g for small k* . The failure of the ansatz for large 
K* reflects the disappearance of the nonlinear effects in 
this range. 

Inserting qcs (as obtained from the wall model) into 
Eq. ^ provides the interaction potential, exhibiting a 
weak dependence on the screening length: 

,a,.).^iic„..(05l.'(§),,i(„-,»).(2) 

Here Xb — (3 e^ j {An €2^0) ~ 0.7 nm is the Bjerrum length 
for water. As discussed before becomes a constant 
of the order of 1 for large k and the K-dependence of 
U is contained only in the wall term U cx ln^[cr*K*^^]. 
The comparison with the predictions of the linear the- 
ory, U oc (cr*K*~^)^, shows that the nonlinear PB the- 
ory yields a drastically changed dependence on both the 
charge density and the screening length. 

Comparison with experiment. There are a number of 
publications concerning the aggregation and compression 
behavior of colloidal mononolayers. However, direct or 
indirect measurements of pair potentials are rather lim- 
ited. Reference [ll| reports results for two batches of 
polystyrene (PS) particles: (a) charged sulfate groups, 
R = 0.55 fim, Uc = 12.5 /zC/cm^ and (b) charged car- 
boxyl groups, R ~ 0.5 /im, Uc — 2.8 /iC/cm^, both inves- 
tigated at an interface between air and ultrapure water 
(k* 1). The tail of the repulsive potential (obtained by 
inverting pair correlation functions) was fitted to a dipole 
form (see Tab. IJ). Comparison with the present renor- 
malized theory (Eq.[2|) requires knowledge of 9. Two dif- 
ferent visual methods applied to sulfonated PS particles 
at the air-water interface 0, [l2| yield quite different re- 
sults which also affects the theoretical result (see Tab. HI. 
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The comparison between experimental and theoretical 
values reveals that for the air-water interface the renor- 
malized charges on the water side seem to be sufficient 
to explain the observed repulsions. In this case charge 
renormalization is essential because the straightforward 
application of the linear theory (Eq. ^) with the bare 
charge gives f3U ~ 10'' x {R/d)^, which is orders of mag- 
nitude off. For the oil-water interface, Ref. reports 
tweezer measurement data for the effective pair potential 
(PS spheres with sulfate groups, R = 1.35 ^m, Cc = 8.9 
liC/cm^) for two electrolyte concentrations {k* sa 2 and 
130). Both data sets could be fitted to one and the 
same pair potential. Although the uncertainty in the 
contact angle translates into a considerable spread of the 
theoretical predictions, the renormalized theory yields a 
potential which is too small by at least a factor of 20. 
Therefore, the experimental results of Ref. @ point to 
still another source of repulsion between the colloids. In 
Ref. Q this other source was argued to be colloidal sur- 
face charges on the oil side, inferred only from the lack 
of a strong K-dependence in the repulsion as predicted 
by the linear theory. This argument is insufficient since 
the renormalized interaction weakens the K-dependence 
considerably; the hypothesis of possible extra charges on 
the oil side is rather supported by the insufficient mag- 
nitude of the renormalized potential. Asymptotically the 
charges on the oil side together with their image charges 
in the water create a net dipole in the nonpolar phase 
Poii ~ 2(jfoii^, where h is their average distance from the 
bottom of the particle. Due to the high dielectric con- 
stant of water, pon will be rather independent of the elec- 
trolyte concentration [l^. The total effective dipole mo- 
ment of the colloid is then given by Pon + Pwatcr where 

J'watcr 

= (2ei/e2)(j'K ^g{a* ,K* ,6) is the dipole moment 
caused by the charges on the water side. The asymp- 
totic interaction between two colloids is then dominated 
by true dipole-dipole interactions given by 

U{d) ~ (Poil + 2poilPwatcr)4T ' 

oTreoei d"^ 

which act in addition to the interaction given in Eq. ([T]). 
The results in Tab. [J suggest that Poii is at least -x/^O 
times larger than Pwater (for pure water). Even then, a 
certain electrolyte concentration dependence of the in- 
teraction potential can be expected through the ensuing 
cross term cx PoiiPwater(K*) which has not been discussed 
in Ref. 0]. We note that recent, more extensive tweezer 
measurements at an oil-water interface show indeed a 
marked dependence on the electrolyte concentration [3] ■ 
In summary, within Poisson-Boltzmann theory we have 
discussed the electrostatic interaction of charged spher- 
ical colloids trapped at an interface between a nonpo- 
lar medium and water. For charges on the water side 
only, we have found a strong renormalization of the ef- 
fective repulsion U , changing the dependence on the sur- 
face charge density Oc and the screening length from 



(7* K* pU/iR/df X 10^ e 

exp. theory (Ref. [5]. .. [12]) 

air /water - Ref. f 11] 
3900 1 8.06 1.8... 4.6 30°... 80° 
800 1 2.16 1.1 ...2.8 30°... 80° 

oil/ water - Ref. [6] 
6800 2 220 1.2 .. . 10 75°. . . 124° 
6800 135 220 0.3 . . . 2.3 75°... 124° 



TABLE I: Comparison between available experimental data 
and Eq. for the amplitude of the interaction potential. For 
simplicity here gun = 1- 



U (X cr^K"^ (linear theory) to [/ cx ln^[e/3/(eoe2)o'c'«^^]. 
For very large charge densities, there is a possibility of a 
near independence of the effective interactions on the salt 
concentration. Geometric effects induced by the shape 
of the colloid are not expected to alter this result signifi- 
cantly as long as is smaller than the linear size of the 
colloid. For colloids at an air-water interface, available 
experimental results compare well with the renormalized 
theory, while for colloids at an oil-water interface the 
renormalized theory underestimates the observed effec- 
tive potential, pointing to an additional source of repul- 
sion such as possible residual charges on the oil side. 
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